-1- 



IMPROVEMENTS IN AND RELATING TO ANALYSIS OF DNA SAMPLES 

This invention concerns improvements in and relating to analysis of DNA 
samples, particularly, but not exclusively, in relation to analysis of DNA samples formed 
of only a few cells. 

A variety of situations, including forensic investigation, make it desirable to be 
able to obtain information about DNA in a sample and express how reliable that 
information is. There are particular problems in analysing small DNA samples and as such 
present techniques tend to involve a substantial level of amplification for the DNA 
followed' by examination of the results by an expat in interpreting such investigations. 
Such exarninations are extremely complex. The expert generally deploys his knowledge to 
determine which of the individual results can be relied upon and which should be 
discounted when reaching the overall result. This approach by its very nature is subjective 
and only makes use of some of the actual individual results obtained when making a 
determination. Additionally the need for an expert analysis introduces a restriction on who 
can perform the review, and hence on the cost and time taken to perform the review. 

The present invention has amongst its aims to provide a technique in which all 
information in an analysis result for a DNA sample is compared with results for reference 
samples with a view to determining a probability of a match between the test sample and 
each of me reference samples by calculating likelihood ratios'. The present invention has 
amongst its aims to provide a technique in which all of the information obtained from the 
analysis of the DNA sample is used in the determination of the results. The present 
invention has amongst its aims to provide a technique in which the impact of potential 
spurious results can be quantified. The present invention has amongst its aims a technique 
for validating approximations which may be made in the analysis of DNA sample. The 
present invention has amongst its aims the provision of a system which can be operated 
successfully by competent but non-expert persons. The present invention has amongst its 
aims the provision of providing a technique whereby complex samples that comprise 
mixtures may be analysed. 
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According to a first aspect of the invention we provide a method of comparing 
one or more reference samples of DNA with at least part of a test sample of DNA, the 
method including:- 

the detennination of the identity of the alleles present at a locus for the DNA in 
the test sample, the determination defining an individual test result, the determination 
being performed for a plurality of loci to give a plurality of individual test results, 

the consideration of one of the plurality of individual test results against the 
individual reference result of one of the reference samples for the respective loci, the 
consideration involving an expression of the probability that the individual reference result 
for that locus could lead by various possible routes to the individual test result for that 
locus, the possible routes to the individual test result including routes where spurious 
information contributes to the individual test result; 

the consideration being repeated for a plurality of loci, the expressions of 
probability that the individual reference result could lead to the individual test result for 
the plurality of loci being combined to give an expression of the probability that the test 
sample matches the reference sample by calculating a likelihood ratio . 

The reference samples may be from known individuals and/or associated with 
other known factors, such as locations, items or events. Each reference sample is 
preferably represented by one or more individual stored results. Each reference sample is 
preferably represented by 8 or more individual stored results. Each reference sample is 
preferably represented by individual stored results which provide the allele identity or 
allele identities for a given locus. 

The reference samples may be stored in a database. The database may be updated 
periodically. New reference samples may be added to the database. Additional individual 
results may be added to existing reference samples. 

The test sample may be from a known individual and/or be associated with one or 
more other known factors, such as a location, item or event the sample was recovered 
from. The test sample may be from one or more sources. One or more of the sources of 
the sample may be known or predicted. 

The identity of the alleles at one or more of the loci of the reference sample and/or 
test sample may be determined by short tandem repeat based investigation. 



Whilst the technique is applicable to all loci, the loci for which allele identity is 
determined may particularly be selected to include one or more of HUMVWFA3 1, 
HUMTH01 3 D21S11, D18S51, HUMFIBRA, D8S1179, HUMAMGXA, HUMAMGY, 
D3S1358, HUMVWA, D16S539, D2S1338, Arnelogenin, D8S1 179, D21S1 1, D1SS51, 
D19S433, HUMTH01 , HUMFBRA/FGA. The loci selected may particularly be each of 
D3S1358, HUMVWA, D16S539, D2S1338, Amelogenin, D8S1 179, D21S1 1, D18S51, 
D19S433, HUMTH01, HUMFIBRA/FGA, 

An individual test result is preferably the allele or alleles detected for that given 
locus and/or the apparent alleles detected when that locus is considered. The individual 
test result may comprise 0 to 50 alleles, but more usually comprises 0 to 32 alleles. The 
individual test result may thus include the homozygous or heterozygous alleles of the test 
sample's source, homozygous and/or heterozygous alleles from DNA contamination of the 
test sample and/or stutters of these and/or other amplification artifacts. 

The consideration is preferably provided for between 1 and 25 loci and more 
preferably between 1 and 16 locL Preferably an equivalent consideration process is 
provided for each locus. 

The consideration preferably involves the determination of a likelihood ratio. The 
likelihood ratio preferably accounts for the probability of the individual sample result 
arising from the individual reference result against the probability of the individual sample 
result arising from other than tide individual reference result 

The consideration may involve the probability of the individual test result arising 
given that individual reference result, including through spurious information occurrence, 
divided by the probability of the individual test result arising from other than the 
individual reference result in any way, including through spurious information occurrence. 
Preferably the consideration may involve the probability of the individual test result 
arising given that individual reference result, including through spurious information 
occurrence, divided by the product of the probability of the individual test result arising 
from other than the individual reference result, including through spurious information 
occurrence, and the frequency of that individual stored result in a population. More 
preferably the consideration may involve the probability of the individual test result arising 
that individual reference result, including through spurious information occurrence, for 
each individual test result, divided by the product of the probability of the individual test 
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result arising from other than the individual reference result, including through spurious 
information occurrence, and the frequency of that individual reference result in a 
population, for each individual test result 

Complex likelihood ratios may be formulated in order to evaluated a mixture. In 
such a case, for a blown and unknown contributor scenario, the likelihood ratio may be the 
probability of the individual test result arising from an individual stored result, and other 
than the individual stored result divided by the probability of the individual test result 
arising from other than the individual stored result and from other than the individual 
stored result. 

The consideration may incorporate an assessment of spurious alleles (either 
stutters of contaminants or other artefacts), that are factored into the probability 
calculations. In addition, the probability of observation of alleles may be calculated from 
the frequency of occurrence in relevant populations and used in the consideration. The 
frequency of occurrence may be derived from an Afro-Caribbean, Asian and white 
Caucasian population. The consideration may include an adjustment to probabilities to 
account fox inbreeding. 

The consideration may take into account more than one route involving spurious 
information and / or more than one type of spurious information. 

Where contamination is necessary to lead to the individual test result the 
probability preferably includes a probability term for spurious allele occurrence. 

Where contamination must not occur to lead to the individual test result the 
probability preferably includes a probability term for spurious allele non-occurrence. 

Where stutter is necessary to lead to the individual test result the probability 
preferably includes a probability term for stutter occurrence. 

Where stutter must not occur to lead to the individual test result the probability 
preferably includes a probability term for stutter non-occurrence. 

Where allele dropout is necessary to lead to the individual test result the 
probability preferably includes a probability term for allele dropout occurrence. 

Where allele dropout must not occur to lead to the individual test result the 
probability preferably includes a probability term for allele dropout non-occurrence. 

Where artifact reporting is necessary to lead to the individual test result the 
probability preferably includes a probability term for artifact reporting occurrence. 
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Where artifact reporting must not occur to lead to the individual test result the 
probability preferably includes a probability term for artifect reporting non-occurrence. 

In the following definitions of probability functions, the probability function may 
include, and ideally is a multiple of, the probability of that possible identity occurring in a 
population. 

Reference to a population may include the world population, a representative 
sample there of, an arbitrary selected population, pseudo-random population, database 
content or other population. 

Where the individual test result has two alleles, the individual reference result has 
two alleles and the individual reference result is a match for the individual test result in 
respect of both alleles, then the probability function may involve, with respect to the 
alleles of the individual reference result, one or more of;- a probability term for one or both 
alleles that drop out did not occur; a probability term for one or both alleles that stutter did 
not occur; a probability term for one or both alleles that spurious alleles did not occur, a 
probability term for one or both alleles that other artifacts did not occur. 

Where the individual test result has two alleles, one being one repeat unit less (or 
4 bases less for a tetrameric locus) than the other, the individual reference result has two 
alleles and the individual reference result is a match for the individual test result in respect 
of the higher allele of the individual test result, but not the lower allele, then the 
probability function may involve, with respect to the alleles of the individual reference 
result, one or more of> a probability term for the matching allele that drop out did not 
occur; a probability term for the non-matching allele that drop out did occur; a probability 
term for stutter of the matching allele that stutter did occur, a probability term for one or 
both alleles of the individual test sample that a spurious allele occurs (preferably with a 
term relating to the frequency of occurrence of that allele in a population); a probability 
term for one or both alleles of the individual test sample that an artifact occurs. 

Where the individual test result has two alleles, one being one repeat unit less (4 
bases less for a tetrameric locus) than the other, the individual reference result has two 
alleles and the individual reference result is a match for the individual test result in respect 
of the lower allele, but not the higher allele, then the probability function may involve, 
with respect to the alleles of the individual reference result, one or more of:-a probability 
term for the matching allele that drop out did not occur; a probability term for the non- 



matching allele that drop out did occur, a probability term for stutter of the matching allele 
that stutter did occur; a probability term for one or both alleles of the individual test 
sample that a spurious allele occurs (preferably with a term relating to the frequency of 
occurrence of that allele in a population); a probability term for one or both alleles of the 
individual test sample that an artifact occurs. 

Where the individual test result has two alleles, the individual reference result has 
two alleles and the individual reference result is a match for the individual test result in 
respect of the lower allele of the individual test result, but not the higher allele, then the 
probability function may involve, with respect to the alleles of the individual reference 
result, one or more of> where the non- matching allele is one repeat unit more (or 4 bases 
more for a tetrameric locus) than the higher allele of the individual test result: a probability 
term for the matching allele that drop out did not occur; a probability term for the non- 
matching allele that drop out did occur; a probability term for stutter of the non-matching 
allele that stutter did occur; a probability term for one or both alleles of the individual test 
sample that a spurious allele occurs (preferably with, a term relating to the frequency of 
occurrence of that allele in a population); a probability term for one or both alleles of the 
individual test sample that an artifact occurs: and where the non-matching allele is not one 
repeat unit more (or 4 bases more for a tetrameric locus) than the higher allele of the 
individual test result one or more of: a probability term for the matching allele that drop 
out did not occur; a probability term for the non-matching allele that drop out did occur; a 
probability term for stutter of the non-matching allele that stutter did not occur; a 
probability term for one or both alleles of the individual test sample that a spurious allele 
occurs (preferably with a term relating to the frequency of occurrence of that allele in a 
population); a probability term for one or both alleles of the individual test sample that an 
artifact occurs. 

Where the individual test result has two alleles, the individual reference result has 
two alleles and the individual reference result is not a match for the individual test result in 
respect of the lower or higher allele of the individual test result, then the probability 
function may involve, with respect to the alleles of the individual reference result, one or 
more of - where one of the non- matching allele is one repeat unit more (or 4 bases more 
for a tetrameric locus) than the higher allele of the individual test result a probability term 
for both the non-matching alleles that drop out did not occur; a probability term for stutter 
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of one of the non-matching alleles that stutter did occur; a probability term for one or both 
alleles of the individual test sample that a spurious allele occurs (preferably with a term 
relating to the frequency of occurrence of that allele in a population); a probability term for 
one or both alleles of the individual test sample that an artifact occurs: and where neither 
the non-matching alleles is one repeat unit more (or 4 bases more for a tetrameric locus) 
than the higher allele of the individual test result, one or more of: a probability term for 
both the non-matching alleles that drop out did occur; a probability term for stutter of the 
non-matching alleles that stutter did not occur; a probability term for one or both alleles of 
the individual test sample that a spurious allele occurs (preferably with a term relating to 
the frequency of occurrence of that allele in a population); a probability term for one or 
both alleles of the individual test sample that an artifact occurs. 

Where the individual test result has two alleles, one being one repeat unit less (4 
bases less for a tetrameric locus) than the other, the individual reference result has one 
allele and the individual reference result is a match for the individual test result in respect 
of the higher allele, then the probability function may involve, with respect to the alleles of 
the individual reference result, one or more of:- a probability term for the matching allele 
that drop out did not occur; a probability term for stutter of the matching allele that stutter 
did occur; a probability term for one or both alleles of the individual test sample that a 
spurious allele occurs (preferably with a term relating to the frequency of occurrence of 
that allele in a population); a probability term for one or both alleles of the individual test 
sample that an artifact occurs. 

Where the individual test result has two alleles, one not being one repeat unit less 
(4 bases less for a tetrameric locus) than the other, the individual reference result has one 
allele and the individual reference result is a match for the individual test result in respect 
of the higher allele, then the probability function may involve, with respect to the alleles of 
the individual reference result, one or more of:- a probability term for the matching allele 
that drop out did not occur; a probability term for stutter of the matching allele that stutter 
did not occur; a probability term for one or both alleles of the individual test sample that a 
spurious allele occurs (preferably with a term relating to the frequency of occurrence of 
that allele in a population); a probability term for one or both alleles of the individual test 
sample that an artifact occurs. 



Where the individual test result has two alleles, the individual reference result has 
one allele and the individual refejsnce result is a match for the individual test result in 
respect of the lower allele, then the probability function may involve, with respect to the 
alleles of the individual reference result, one or more of:- a probability term for the" 
matching allele that drop out did not occur; a probability term for stutter of the matching 
allele that stutter did not occur; a probability term for one or both alleles of the individual 
test sample that a spurious allele occurs (preferably with a term relating to the frequency of 
occurrence of that allele in a population); a probability term for one or both alleles of the 
individual test sample that an artifact occurs. 

Where the individual test result has two alleles, the individual reference result has 
one allele and the individual reference result is not a match for the individual test result in 
respect of either allele, then the probability function may involve, with respect to the 
alleles of the individual reference result, one or more of:- where the non-matching allele is 
one repeat unit more (4 bases more for a tetrameric locus) than one of the individual test 
result alleles: a probability term for the non-matching allele that drop out did occur, a 
probability term for stutter of the non-matching allele that stutter did occur, a probability 
term for one or both alleles of the individual test sample that a spurious allele occurs 
(preferably with a term relating to the frequency of occurrence of that allele in a 
population); a probability term for one or both alleles of the individual test sample that an 
artifact occurs: and where the non-matching allele is not one repeat unit more (4 bases 
more for a tetrameric locus) than one of the individual test result alleles one or more of; a 
probability term for the non-matching allele that drop out did occur; a probability term for 
stutter of the non-matching allele thai stutter did not occur; a probability tenn for one or 
both alleles of the individual test sample that a spurious allele occurs (preferably with a 
term relating to the frequency of occurrence of that allele in a population); a probability 
term for one or both alleles of the individual test sample that an artifact occurs . 

Where the individual test result has one allele, the individual reference result has 
two alleles and the lower allele of the individual reference result is a match for the 
individual test result, then the probability function may involve, with respect to the alleles 
of the individual reference result, one or more of:- where the higher allele of the individual 
reference result is one repeat unit more (4 bases more for a tetrameric locus) than the allele 
of the individual test result; a probability term for the matching allele that drop out did not 
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occw> a probability for the non-matching allele that drop out did occur; a probability term 
for stutter of the non-matching allele that stutter did occur; a probability term for one or 
both alleles of the individual test sample that a spurious allele does not occur (preferably 
with a term relating to the frequency of occurrence of that allele in a population); a " 
probability term for one or both alleles of the individual test sample that an artifact does 
not occur: and where the higher allele of the individual reference result is not one repeat 
unit more (4 bases more for a tetrameric locus) than the allele of the individual test result 
one or more of: a probability term for the matching allele that drop out did not occur; a 
probability for the non-matching allele that drop out did occur; a probability term for 
stutter of the non-matching allele that stutter did not occur; a probability term for one or 
both alleles of the individual test sample that a spurious allele does not occur (preferably 
with a term relating to the frequency of occurrence of that allele in a population); a 
probability term for one or both alleles of the individual test sample that an artifact does 
not occur. 

Where the individual test result has one allele, the individual reference result has 
two alleles and neither of the individual reference results is a match for the individual test 
result, then the probability function may involve, with respect to the alleles of the 
individual reference result one or more of> where one of the individual reference result is 
one repeat unit more (4 bases more for a tetrameric locus) than the allele of the individual 
test result; a probability term for one or both the non-matching alleles that drop out did 
occur; a probability for stutter of the non-matching allele which is one repeat unit more (4 
bases more for a tetrameric locus) than the individual test result that stutter did occur; a 
probability term for one or both alleles of the individual test sample that a spurious allele 
does occur (preferably with a term relating to the frequency of occurrence of that allele in a 
population); a probability term for one or both alleles of the individual test sample that an 
artifact does occur: and where neither of the alleles of the individual reference result is one 
repeat unit more (4 bases more for a tetrameric locus) than the allele of the individual test 
result one or more of: a probability term fox both the non-matching alleles that drop out did 
occur; a probability term for stutter of the non-matching allele that stutter did not occur; a 
probability term for one or both alleles of the individual test sample that a spurious allele 
does occur (preferably with a term relating to the frequency of occurrence of that allele in a 
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population); a probability term for one or both alleles of the individual test sample that an 
artifact does occur. 

Where the individual test result has one allele, the individual reference result has 
one allele and the alleles match then the probability function may involve, with respect to 
the alleles of the individual reference result, one or more of:-a probability tenn for the 
matching allele that drop oat did not occur; a probability term for stutter for the matching 
allele that stutter does not occur; a probability term for the allele of the individual test 
sample that a spurious allele does not occur (preferably with a term relating to the 
frequency of occurrence of that allele in a population); a probability tenn for the allele of 
the individual test sample that an artifact does not occur. 

Where the individual test result has one allele, the individual reference result has 
one allele and the alleles do no match then the probability function may involve, with 
respect to the alleles of the individual reference result, one or more of> where the non- 
matching allele of the individual reference result is one repeat unit more (4 bases more for 
a tetrameric locus) than the allele of the individual test result; a probability term for the 
non-matching allele that drop out occurs; a probability tenn for stutter of the non-matching 
allele that stutter occurs; a probability term for the allele of the individual test sample that 
a spurious allele occurs (preferably with a term relating to the frequency of occurrence of 
that allele in a population); a probability tenn for the allele of the individual test sample 
that an artifact does occur. 

The various possible routes for the individual stored result giving the individual 
sample result may include contamination giving one or more alleles in the individual 
sample result not present in the individual stored result 

The various possible routes for the individual stored result giving the individual 
sample result may include stutter giving one or more alleles in the individual sample result 
not present in the individual stored result. 

The various possible routes for the individual stored result giving the individual 
sample result may include amplification of artifacts giving one or more alleles in the 
individual sample result not present in the individual stored result. 

The various possible routes for the individual stored result giving the individual 
sample result may include allele drop out giving one or more alleles missing in the 
individual sample result present in the individual stored result 
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The probability function may include a probability that contamination may occur. 
The probability that contamination may occur may be determined by one or more control 
determinations. The control determinations may be made in parallel with the 
determination of the identity of the alleles of the test sample. The control determinations 
may be made separately, for instance as a reference investigation used subsequently in two 
or more test sample determinations. The probability that contamination may occur may be 
provided for by theoretical predictions. 

Hie probability function may include a probability that stutter may occur. The 
probability that stutter may occur may be determined by one or more control 
determinations. The control detenninations may be made in parallel with the 
determination of the identity of the alleles of the test sample. The control determinations 
may be made separately, for instance as a reference investigation used subsequently in two 
or more test sample determinations. The probability that stutter may occur may be 
provided for by theoretical predictions. 

The probability function may include a probability that allele dropout may occur. 
The probability that allele dropout may occur may be determined by one or more control 
determinations. The control detenninations may be made in parallel with the 
determination of the identity of the alleles of the test sample, The control detenninations 
may be made separately, for instance as a reference investigation used subsequently in two 
or more test sample determinations. The probability that allele dropout may occur may be 
provided for by theoretical predictions. 

The probability function may include a probability that artifact reporting may 
occur. The probability that artifact reporting may occur may be determined by one or 
more control determinations. The control determinations may be made in parallel with the 
determination of the identity of the alleles of the test sample. The control determinations 
may be made separately, for instance as a reference investigation used subsequently in two 
or more test sample detenninations. The probability that artifact reporting may occur may 
be provided for by theoretical predictions. 

The spurious information may be due to contamination effects, allele dropout 
effects, locus dropout effects, stutter effects, artifact effects or other causes. 

The contribution of the spurious information may lead to an allele being present 
which is not part of the DNA test sample, the absence of alleles which should be present 
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from the DNA test sample, the presence of apparent alleles in positions one repeat unit (4 
bases lower for a tetrameric locus) than the alleles in the DNA test sample. 

Preferably the consideration is applied to a plurality of loci 7 ideally all loci for 
which individual stored results and/or individual test results exist. 

The combination of probabilities produced by the respective considerations is 
preferably obtained by multiplying the probabilities together. 

Two or more different determinations of the identities of the alleles in the test 
sample may be performed. The method may be applied to each set of individual test 
results thereby obtained. The expression of a likelihood ratio for respective sets of 
individual test results may be considered against one another and/or combined. 

The expression of a likelihood ratio and/or a combined expression of a likelihood 
ratio that a given reference sample and test sample match may be generated for a plurality, 
ideally all, of the reference samples available. The reference samples may be ranked in 
order of the likelihood ratios of a match with the test sample, ideally descending order. 

According to a second aspect of the invention we provide a method of indicating a 
likelihood ratio that evaluates that at least a part of a DNA test sample arose from a known 
source, the method involving:- 

one or more determinations of the identity of the alleles present at a locus for the 
DNA in the test sample, each determination defining an individual test result; 

the determination of at least some of the theoretical allele identities which could 
have produced a given individual test result, these identities forming the individual 
reference results; 

the determination of the identity of the alleles present at the locus for the DNA 
from the known source; 

one of the theoretical allele identities being the identity determined for that locus 
for the known source; 

the provision of a probability function for each individual reference result 
considered which is representative of at least some of the various possible routes by which 
that given individual reference result may lead to the given individual test result, that 
probability function further being representative of the likelihood of that individual 
reference result's occurrence and the possible routes to the individual test result which 
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includes routes where spurious information contributes, this probability function forming 
defining the theoretical probability functions; 

the theoretical probability functions for different individual reference results 
being combined to give an indication of the various ways in which the given individual test 
result could be reached, this combination forming the combined theoretical probability 
function; 

the provision of a probability function for the individual reference result matching 
the known source's identity, which is representative of the manner in which that individual 
reference result leads to the individual test result, this forming the known source's 
theoretical function* 

the known source's theoretical function and combined theoretical function being 
considered together to calculate the likelihood ratio. 

The second aspect of the invention may include features, options or possibilities 
set out elsewhere in this document. 

At least part of a DNA sample may refer to one source of a multi-source or mixed 
sample. The method may indicate calculation of a likelihood ratio relating to one or more 
sources of a defined nature, for instance the likelihood ratio may evaluate the proposition 
of two defined contributors to the sample. 

The known source may refer to a known individual and/or be associated with one 
or more other known factors, such as a location, item or event the sample was recovered 
from. 

The identity of the alleles at one or more of the loci of the test sample may be 
determined by short tandem repeat based investigation. 

An individual sample result is preferably the allele or alleles detected for that 
given locus and/or the apparent alleles detected when that locus is considered. The 
individual sample result may comprise 0 to 50 alleles but more usually comprises 0 to 32 
alleles. The individual sample result may thus include lie homo2ygous or heterozygous 
alleles of Hhe test sample's source, homozygous and/or heterozygous alleles from DNA 
contamination of the test sample and/or stutters of these and/or other amplification 
artifacts- 
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The consideration is preferably provided for between 1 and 25 loci and more 
preferably between 1 and 16 loci. Preferably an equivalent consideration process is 
provided for each loci. 

Whilst the technique is applicable to all loci, the loci for which allele identity is 
determined may particularly be selected to include one or more of HUMVWFA3 1, 
HUMTH01,D21Sll ? D18S51, HUMFffiRA, D8SU79, HUMAMGXA, HUMAMGY, 
D3S1358, HUMVWA, D16S539, D2S1338, Amelogenin, D8S1179, D21S11, D18S51. 
D19S433, HUMTH01, HUMFIBRA/FGA. The loci selected may particularly be each of 
D3S1358, HUMVWA, D16S539 7 D2S1338, Amelogenin, D8S1179, D21S11, D18S51, 
D19S433, HUMTH01, HUMFIBRA/FGA. 

The theoretical identities may be determined from the alleles indicated in the 
individual test result, All possible theoretical identities may be determined, but more 
preferably those theoretical identities which could reasonably lead to the individual test 
result are determined Those theoretical identities defined as reasonable may be all 
identities where an allele in the test sample is in common with the reference sample. The 
determination may involve providing theoretic identities corresponding to each 
permutation of two alleles, where at least one of those alleles matches an allele in the 
individual test result. 

The provision of a theoretical probability function may involve the probability of 
getting that individual test result in any way, including through spurious information 
occurrence. Preferably the provision of a probability function may involve the probability 
of getting that individual test result in any way, including through spurious information 
occurrence, and the frequency of that given theoretical identity in a population. More 
preferably the provision of a probability function may involve the probability of getting 
that individual test result in any way, including through spurious information occurrence, 
and the frequency of that theoretical identity in a population, for each individual test result. 

The theoretical probability function for each individual reference result theoretical 
identity is preferably defined in part by a probability for that individual reference results 
identity occurrence in a population. The theoretical probability function for each 
individual reference result is preferably defined in part by a probability for the various 
occurrences which would result in that individual reference result giving the individual test 
result 
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Theoretical probability functions may be provided to account for each of the 
individual test results determined for a locus in the aforementioned manner. Preferably the 
theoretical probability functions for each individual test result given an individual 
reference result are combined, ideally before the theoretical probability function s for 
different individual reference results are combined. Preferably the theoretical probability 
functions for different individual test results are combined by multiplication. Preferably 
the theoretical probability functions for different individual reference results are combined 
by addition. 

Where contamination is necessary to lead to the individual test result the 
probability preferably includes a probability term for spurious allele occurrence. 

Where contamination must not occur to lead to the individual test result the 
probability preferably includes a probability term for spurious allele non-occurrence. 

Where stutter is necessary to lead to the individual test result the probability 
preferably includes a probability term for stutter occurrence. 

Where stutter must not occur to lead to the individual test result the probability 
preferably includes a probability term for stutter non-occurrence. 

Where allele dropout is necessary to lead to the individual test result the 
probability preferably includes a probability term for allele dropout occurrence. 

Where allele dropout must not occur to lead to the individual test result the 
probability preferably includes a probability term for allele dropout non-occurrence. 

Where artifact reporting is necessary to lead to the individual test result the 
probability preferably includes a probability term for artifact reporting occurrence. 

Where artifact reporting must not occur to lead to the individual test result the 
probability preferably includes a probability term for artifact reporting non-occurrence. 

In the following definitions of probability functions, the probability function may 
include, and ideally is a multiple of, the probability of that possible identity occurring in a 
population. 

Reference to a population may include the world population, a representative 
sample there of, an arbitrary selected population* pseudo-random population, database 
content or other population. 

Where the individual test result has two alleles, the individual reference result has 
two alleles and the individual reference result is a match for the individual test result in 
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respect of both alleles ; then the probability function may involve, with respect to the 
alleles of the individual reference result, one or more of:- a probability term for one or both 
alleles that drop out did not occur; a probability term for one or both alleles that stutter did 
not occur, a probability term for one or both alleles that spurious alleles did not occur, a 
probability term for one or both alleles that other artifacts did not occur. 

Where the individual test result has two alleles, one being one repeat unit less (or 
4 bases less for a tetrameric locus) than the other, the individual reference result has two 
alleles and the individual reference result is a match for die individual test result in respect 
of the higher allele of the individual test result, but not the lower allele, then the 
probability function may involve, with respect to the alleles of the individual reference 
result, one or more of:- a probability term for the matching allele that drop out did not 
occur; a probability term for the non-matching allele that drop out did occur; a probability 
term for stutter of the matching allele that stutter did occur; a probability term for one or 
both alleles of the individual test sample that a spurious allele occurs (preferably with a 
term relating to the frequency of occurrence of that allele in a population); a probability 
term for one or both alleles of die individual test sample that an artifact occurs. 

Where the individual test result has two alleles, one being one repeat unit less (4 
bases less for a tetrameric locus) than the other, the individual reference result has two 
alleles and the individual reference result is a match for the individual test result in respect 
of tie lower allele, but not the higher allele, then the probability function may involve, 
with respect to the alleles of the individual reference result, one or more of> a probability 
term for the matching allele that drop out did not occur; a probability term for the non- 
matching allele that drop out did occur; a probability term for stutter of the matching allele 
that stutter did occur, a probability term for one or both alleles of the individual test 
sample that a spurious allele occurs (preferably with a term relating to the frequency of 
occurrence of that allele in a population); a probability term for one or both alleles of the 
individual test sample that an artifact occurs. 

Where the individual test result has two alleles, the individual reference result has 
two alleles and the individual reference result is a match for the individual test result in 
respect of the lower allele of the individual test result, but not the higher allele, then the 
probability function may involve, with respect to the alleles of the individual reference 
result, one or more of:- where the non- matching allele is one repeat unit more (or 4 bases 
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more for a tetrameric locus) than the higher allele of the individual test result, a probability 
term for the matching allele that drop out did not occur; a probability term for the non- 
matching allele that drop out did occur, a probability term for stutter of the non-matching 
allele that stutter did occur; a probability term for one or both alleles of the individual test 
sample that a spurious allele occurs (preferably with a term relating to the frequency of 
occurrence of that allele in a population); a probability term for one or both alleles of the 
individual test sample that an artifect occurs: and where the non-matching allele is not one 
repeat unit more (or 4 bases more for a tetrameric locus) than the higher allele of the 
individual test result on or more of: a probability term for the matching allele that drop out 
did not occur; a probability term for the non-matching allele that drop out did occur; a 
probability term for stutter of the non-matching allele that stutter did not occur, a 
probability terra for one or both alleles of the individual test sample that a spurious allele 
occurs (preferably with a term relating to the frequency of occurrence of that allele in a 
population)- a probability term for one or both alleles of the individual test sample that an 
artifact occurs. 

Where the individual test result has two alleles, the individual reference result has 
two alleles and the individual reference result is not a match for the individual test result in 
respect of the lower or higher allele of the individual test result, then the probability 
function may involve, with respect to the alleles of the individual reference result, one or 
more of:- where one of the non- matching allele is one repeat unit more (or 4 bases more 
for a tetrameric locus) than the higher allele of the individual test result, a probability term 
for both the non-matching alleles that drop out did not occur; a probability term for stutter 
of one of the non-matching alleles that stutter did occur, a probability term for one or both 
alleles of the individual test sample that a spurious allele occurs (preferably with a term 
relating to the frequency of occurrence of that allele in a population); a probability tenn for 
one or both alleles of the individual test sample that an artifact occurs; and where neither 
the non-matchiDg alleles is one repeat unit more (or 4 bases more for a tetrameric locus) 
than the higher allele of the individual test result one or more of: a probability term for 
both the non-matching alleles that drop out did occur, a probability term for stutter of the 
non-matching alleles that stutter did not occur; a probability term for one or both alleles of 
the individual test sample that a spurious allele occurs (preferably with a term relating to 
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the frequency of occurence of that allele in a population); a probability term for one or 
both alleles of the individual test sample that an artifact occurs. 

Where the individual test result has two alleles, one being one repeat unit less (4 
bases less for a tetrameric locus) than the other, die individual reference result has'one 
allele and the individual reference result is a match for the individual test result in respect 
of the higher allele, then the probability function may involve, with respect to the alleles of 
the individual reference result, one or more of:- a probability term for the matching allele 
that drop out did not occur; a probability term for stutter of the matching allele that stutter 
did occur; a probability term for one or both alleles of the individual test sample that a 
spurious allele occurs (preferably with a term relating to the frequency of occurrence of 
that allele in a population); a probability term for one or both alleles of the individual test 
sample that an artifact occurs. 

Where the individual test result has two alleles, one not being one repeat unit less 
(4 bases less for a tetrameric locus) than the other, the individual reference result has one 
allele and the individual reference result is a match for the individual test result in respect 
of the higher allele, then the probability function may involve, with respect to the alleles of 
the individual reference result, one or more of> a probability term for the matching allele 
that drop out did not occur; a probability term for stutter of the matching allele that stutter 
did not occur; a probability term for one or both alleles of the individual test sample that a 
spurious allele occurs (preferably with a term relating to the frequency of occurrence of 
that allele in a population); a probability term for one or both alleles of the individual test 
sample that an artifact occurs. 

Where the individual test result has two alleles, the individual reference result has 
one allele and the individual reference result is a match for the individual test result in 
respect of the lower allele, then the probability function may involve, with respect to the 
alleles of the individual reference result, one or more of:- a probability term for the 
matching allele that drop out did not occur; a probability term for stutter of the matching 
allele that stutter did not occur; a probability term for one or both alleles of the individual 
test sample that a spurious allele occurs (preferably with a term relating to the frequency of 
occurrence of that allele in a population); a probability term for one or both alleles of the 
individual test sample that an artifact occurs. 
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Where die individual test result has two alleles, the individual reference result has 
one allele and the individual reference result is not a match for the individual test result in 
respect of either allele, then the probability function may involve, with respect to the 
alleles of the individual reference result, one or more of> where the non-matching allele is 
one repeat unit more (4 bases more for a tetrameric locus) than one of the individual test 
result alleles, a probability term for the non-matching allele that drop out did occur; a 
probability term for stutter of the non-matching allele that stutter did occur; a probability 
term for one or both alleles of the individual test sample that a spurious allele occurs 
(preferably with a term relating to die frequency of occurrence of that allele in a 
population); a probability term for one or both alleles of the individual test sample that an 
artifact occurs: and where the non-matching allele is not one repeat unit more (4 bases 
more for a tetrameric locus) than one of the individual test result alleles, on or more of; a 
probability term for the non-matching allele that drop out did occur; a probability term for 
stutter of the non-matching allele that stutter did not occur, a probability term for one or 
both alleles of the individual test sample that a spurious allele occurs (preferably with a 
term relating to the frequency of occurrence of that allele in a population); a probability 
term for one or both alleles of the individual test sample that an artifoct occurs . 

Where the individual test result has one allele, the individual reference result has 
two alleles and the lower allele of the individual reference result is a match for the 
individual test result, then the probability function may involve, with respect to the alleles 
of the individual reference result, one or more of> where the higher allele of the individual 
reference result is one repeat unit more (4 bases more for a tetrameric locus) than the allele 
of the individual test result: a probability term for the matching allele that drop out did not 
occur, a probability for the non-matching allele that drop out did occur; a probability term 
for stutter of the non-matching allele that stutter did occur; a probability term for one or 
both alleles of the individual test sample that a spurious allele does not occur (preferably 
with a term relating to the frequency of occurrence of that allele in a population); a 
probability term for one or both alleles of the individual test sample that an artifact does 
not occur: and where the higher allele of the individual reference result is not one repeat 
unit more (4 bases more for a tetrameric locus) than the allele of the individual test result, 
one or more of: a probability term for the matching allele that drop out did not occur, a 
probability for the non-matching allele that drop out did occur; a probability term for 
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stutter of the non-matching allele that stutter did not occur; a probability term for one or 
botii alleles of the individual test sample that a spurious allele does not occur (preferably 
with a term relating to the frequency of occurrence of that allele in a population); a 
probability term for one or both alleles of the individual test sample that an artifacr does 
not occur. 

Where the individual test result has one allele, the individual reference result has 
two alleles and neither of the individual reference results is a match for the individual test 
result, then the probability function may involve, with respect to the alleles of die 
individual reference result, one or more of> where one of the individual reference result is 
one repeat unit more (4 bases more for a tetrameric locus) than the allele of the individual 
test result: a probability term for one or both the non-matching alleles that drop out did 
occur; a probability for stutter of the non-matching allele which is one repeat unit more (4 
bases more for a tetrameric locus) than the individual test result that stutter did occur; a 
probability term for one or both alleles of the individual test sample that a spurious allele 
does occur (preferably with a term relating to the frequency of occurrence of that allele in a 
population); a probability term for one or both alleles of the individual test sample that an 
artifact does occur: and where neither of the alleles of the individual reference result is one 
repeat unit more (4 bases more for a tetrameric locus) than the allele of the individual test 
result, one or more of: a probability term for both the non-matching alleles that drop out 
did occur; a probability term for stutter of the non-matching allele that stutter did not 
occur, a probability term for one or both alleles of the individual test sample that a 
spurious allele does occur (preferably with a term relating to the frequency of occurrence 
of that allele in a population); a probability terra for one or both alleles of the individual 
test sample that an artifact does occur. 

Where the individual test result has one allele, the individual reference result has 
one allele and the alleles match then the probability function may involve, with respect to 
the alleles of the individual reference result, one or more of> a probability term for the 
matching allele that drop out did not occur; a probability term for stutter for the matching 
allele that stutter does not occur; a probability term for the allele of the individual test 
sample that a spurious allele does not occur (preferably with a term relating to the 
frequency of occurrence of that allele in a population); a probability term for the allele of 
the individual test sample that an artifact does not occur. 
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Where the individual test result has one allele, the individual reference result has 
one allele and the alleles do no match then the probability function may involve, with 
respect to the alleles of the individual reference result, one or more of:- where the non- 
matching allele of the individual reference result is one repeat unit more (4 bases more for 
a tetrameric locus) than the allele of the individual test result: a probability term for the 
non-matching allele that drop out occurs; aprobability term for stutter of the non-matching 
allele that stutter occurs; a probability term for the allele of the individual test sample that 
a spurious allele occurs (preferably with a term relating to the frequency of occurrence of 
that allele in a population); a probability terra for the allele of the individual test sample 
that an artifact does occur. 

The various possible routes for the individual reference result giving the 
individual test result may include contamination giving one or more alleles in the 
individual test result not present in the individual reference result 

The various possible routes for the individual reference result giving the 
individual test result may include stutter giving one or more alleles in the individual test 
result not present in the individual reference result 

The various possible routes for Hie individual reference result giving the 
individual test result may include amplification of artifacts giving one or more alleles in 
the individual test result not present in the individual reference result. 

The various possible routes for the individual reference result giving the 
individual test result may include allele drop out giving one or more alleles missing in the 
individual test result present in the individual reference result. 

The probability function may include a probability that contamination may occur. 
The probability that contamination may occur may be determined by one or more control 
determinations. The control determinations may be made in parallel with the 
determination of the identity of the alleles of the test sample. The control determinations 
may be made separately, for instance as a reference investigation used subsequently in two 
or more test sample determinations. The probability that contamination may occur may be 
provided for by theoretical predictions. 

The probability function may include a probability that stutter may occur. The 
probability that stutter may occur may be determined by one or more control 
determinations. The control determinations may be made in parallel with the 
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detennination of the identity of the alleles of the test sample. The control determinations 
may be made separately, for instance as a reference investigation used subsequently in two 
or more test sample determinations. The probability that stutter may occur may be 
provided for by theoretical predictions. 

The probability function may include a probability that allele dropout may occur. 
The probability that allele dropout may occur may be determined by one or more control 
determinations. The control determinations may be made in parallel with the 
determination of the identity of the alleles of the test sample. The control determinations 
may be made separately, for instance as a reference investigation used subsequently in two 
or more test sample determinations, The probability that allele dropout may occur may be 
provided for by theoretical predictions. 

The probability function may include a probability that artifact reporting may 
occur. The probability that artifact reporting may occur may be determined by one or 
more control determinations. The control determinations may be made in parallel with the 
determination of the identity of the alleles of the test sample. The control determinations 
may be made separately, for instance as a reference investigation used subsequently in two 
or more test sample determinations. The probability that artifact reporting may occur may 
be provided for by theoretical predictions. 

The spurious information may be due to contamination effects, allele dropout 
effects, locus dropout effects, stutter effects, artifact effects or other causes. 

The contribution of the spurious information may lead to an allele being present 
which is not part of the DNA test sample, the absence of alleles which should be present 
from the DNA test sample, the presence of apparent alleles in positions one repeat unit less 
(or 4 bases less for a tetrameric locus) than the alleles in the DNA test sample. 

The theoretical probability functions may be combined to give the overall 
combined theoretical probability function by summing the theoretical probability functions 
together. 

The provision of the probability function results matching the known source's 
identity may involve the probability of getting that individual test result given that 
individual reference result, including through spurious information occurrence. Preferably 
the provision of the probability function for the individual reference may involve the 
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probability of getting that individual test result given that individual reference result, 
including through spurious information occurrence, for each individual test result 

The known source's theoretical function and combined theoretical function may 
be combined as a ratio, preferably as a likelihood ratio. The likelihood ratio preferably 
accounts for the probability that a given individual reference result / theoretical identity 
leads to the individual test result against the probability that the individual test result was 
lead to in another way. The likelihood ratio may be the known source's theoretical 
function divided by the combined theoretical function. 

Preferably the method is repeated for a plurality of loci, ideally all loci for which 
individual test results exist. The likelihood ratio obtained for each loci may be multiplied 
together to give a combined loci likelihood ratio. 

Two or more different determinations of the identities of the alleles in the test 
sample may be performed. The method may be applied to each set of individual test 
results thereby obtained. The expression of the likelihood ratio for respective sets of 
individual test results may be considered against one another and/or combined. 

According to a third aspect of the invention we piovide a method of investigating 
the acceptable values for one or more variables relating to a DNA sample analysis, the 
method involving:- 

the consideration of one or more probability functions used in a method of 
indicating the likelihood ratio that at least part of a DNA sample arose from a known 
source, at least one of the probability functions being defined by an approximating 
function and a scaling factor, the scaling factor including at least one of the variables as a 
term; 

the value for one or more of the variables in the scaling fector being assigned a 
plurality of different values and the value of the scaling factor being considered at each of 
those different values; 

the value or values for the one or more variables being deemed acceptable when 
the value of the scaling fector is within a predetermined or acceptable range. 
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The third aspect of the invention may include any of the features, options or 
possibilities set out elsewhere in this document, including the first and/or second aspects of 
the invention. 

The DNA sample analysis is preferably a consideration of the likelihood ratio that 
a sample arose from one or more scenarios compared with the sample arising from the 
other possible scenarios. 

The probability functions may be particularly provided according to the 
definitions of the first and/or second aspects of the invention. The probability functions 
preferably take into account the probability of spurious information potentially 
contributing to the results obtained upon analysis of the sample. 

The approximating function may provide an indication of the probability within 
certain acceptable ranges for potential variables of the analysis process. The 
approximating function may be an accurate assumption of the probability within these 
acceptable ranges. The approximating function may be an inaccurate assumption of the 
probability outside these acceptable ranges. 

The scaling factor may account for one or more variables in the analysis. The one 
or more variables may be sources of error, The one or more variables may be probabilities 
of spurious information contributing to the results of the analysis of the sample. The 
spurious information sources may be one or more of allele/locus dropout, stutter, 
contamination or artifact reporting. 

The variable value is preferably between 0 and 1 inclusive. The variable value 
may be assigned values at increments of 0.1 or less for the investigation. Preferably each 
combination of the incremental values is considered for the variables which contribute to a 
given scaling factor. A scaling factor may involve one, two, three, four or more variables 
depending upon the scaling factor. Preferably values for the scaling factor are determined 
for all possible combinations of the variable's values. 

The acceptable range is preferably a range in which the scaling factor has a 
minimal effect on the probability function if it is included compared with if it is excluded. 
The range for the scaling factor may be between 0.9 and 1.1 in some cases. In other cases 
the scaling factor may be between 0.9 and L 
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According to a fourth aspect of the invention we provide a method of indicating a 
likelihood ratio that at least a part of a DNA sample arose from a known source or sources, 
the method involving;* 

the determination of the identity of the alleles present at a locus for the DNA in 
the sample, the determination defining an individual test result; 

the determination of the identity of the alleles present at the locus for the known 

source; 

the consideration of a likelihood ratio that the known source leads to the 
individual test result compared with the other possible routes to the individual test result, 
the likelihood ratio being based on ond or more probability functions, at least one of the 
probability functions being defined by an approximating function and a scaling factor, the 
scaling factor including at least one variable term relating to the probability of spurious 
information potentially contributing to the individual result; 

the value of one or more of the variables being determined for the method; 

the determined value of the one oj more variables being considered against an 
acceptable range for that variable and/cjr the value of one or more of the scaling factors 
being considered against its acceptable range given that determined value for that variable, 

one or more of the probability functions defined by a scaling factor including that 
variable being deemed defined by the approximating function where that variable has a 
value within its acceptable range and/or where the scaling factor has a value within its 
predetermined or acceptable range, the so defined one or more probability functions being 

i 

used as the basis for the likelihood ratK>, 

The fourth aspect of the invention may include any of the features, options or 
possibilities set out elsewhere in this document, including the first and/or second and/or 

third aspects of the invention. ' 

i 

The value of one or more of the variables may be determined for the method in a 
different way from the determination of one or more other variables. The determination 
may be carried out for the laboratory where the analysis of the DNA sample will be 
performed. The determination may be performed alongside the analysis. The 

determination may be performed separately form the analysis, including periodic 

1 

determinations and even one off determinations. The determination may be made using 



-26- 

control analyses, including negative and/or positive controls. Experimental deterioration is 
preferred for the contamination value, for instance, The determination may be made 
theoretically. 

The acceptable range is preferably a range in which the scaling factor has a 
minimal effect on the probability function if it is included compared with if it is excluded. 
The range for the scaling factor may be between 0.9 and 1.1 in some cases. In other cases 
the scaling factor may be between 0.9 and 1. 

The probability functions may be defined by the approximating function and 
scaling factor where the variable values for that probability function and/or the scaling 
factor value for that probability function is outside the acceptable ranges. One or more 
probability functions defined by the approximating function and scaling factor may be 
used in combination with one or more probability functions defined by the approximating 
function, but in general all the probability functions will be defined by either the 
approximating function and scaling factor or by the approximating factor alone. 

The invention will now be described, by way of example only, and with reference 
to the accompanying drawings in which :- 

Figure 1 illustrates the consideration of results from DHA 
analysis in a technique not according to the present invention; 

Figure 2 illustrates the consideration of results from DNA 
analysis in a technique according to one embodiment of the present 
invention; 

Figure 3 provides Table 1 which illustrates the calculation of the 
components of the likelihood ratio for an example where three individual 
results show evidence of spurious bands and allele dropout; 

Figure 4 illustrates Table 2 and the calculation of the 
components of the likelihood ratio for an example involving stutter; 

Figure 5 illustrates graphically the testing of the robustness of 
the F designation from example la; 

Figure 6a illustrates graphically the evaluation of scaling factor 
in example 1c; 
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Figure 6b illustrates graphically the evaluation of the scaling 
factor of example lc three dimensionally; 

Figure 7b illustrates graphically the evaluation of the scaling 
function from example 2a; - 

Figure 7b illustrates graphically the evaluation of the scaling 
factor from example 2a three dimensionally; 

Figure 8a illustrates graphically the evaluation of the scaling 
factor from example 2b; 

Figure 8b illustrates graphically the evaluation of the scaling 
factor from example 2b three dimensionally; 

Figure 9 illustrates graphically the evaluation of the scaling 
function from example 3a; and 

Figure 10 represents Table 3, an analysis of p(D) parameters 
derived from experimental observation. 



In any DNA analysis technique the individual results can include information not 
arising from the DNA actually under investigation and/or not accurately represent the 
variations in the DNA in question which is actually being considered. This is particularly 
true for analysis techniques where amplification occurs, such as PCR based analysis. The 
problems become increasingly significant as the size of the initial DNA sample to be 
investigated becomes smaller. By the time investigations on a small number of cells are 
considered, 10 or less for instance, then substantial potential for such issues exists. In 
particular, even the lowest levels of DNA contamination from sample collection, 
laboratory handling or the equipment itself can have a marked effect on the individual 
results obtained. Additionally locus or allele dropout can occur and the amplified products 
may not fully reflect the DNA of the sample. There is also the possibility that stutters in 
the results will be in distinguishable from genuine allele results. 

The outcome of these factors and others is that the potential for analysing small 
DNA samples is presently limited and where analysis is carried out it has to be handled by 
experts. Based on their knowledge and experience such expert individuals are able to 
consider the individual results and discard an individual result from further consideration 
where mere is a question mark over it. This may lead to a substantial number of occasions 
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where the individual result has to be discounted due to a negative control providing a 
response, with consequential cost and time implications. Additionally such processes are 
by their very nature subjective and equivalent sets of individual results may be handled by 
different experts in different ways. The discarding of many individuals from the initial set 
of results from further consideration also means that a high proportion of the individual 
results do not contribute in anyway to the overall findings. 

The techniques of die present invention, whilst particularly useful in relation to 
consideration of small DNA samples, are useful in all DNA analysis techniques based 
around similar considerations. 

The present invention in one embodiment below, aims to provide a system which 
takes in to account all the information obtained on a sample of unknown source and 
expresses a likelihood ratio that it matches with each of various known reference samples, 
rather than indicating a match or non-match against each. The present invention thus aims 
to reduce the skill required to perform an analysis but improve the accuracy and/or level of 
information provided by the analysis. 

The present invention in one embodiment below, aims to provide a system which 
accounts for the other sources of information effectively and in a standardised way. The 
present invention, as embodied below, aims to provide results which are based on all of the 
initial individual results obtained, rather than the results after a screening process. 

In another embodiment, described below, the present invention aims to investigate 
and determine the impact of the other sources of information and particularly their 
probability of occurrence on a process for investigating DNA samples. 

In another embodiment of the present invention, described below, the invention 
aims to provide a system which provides results for a DNA analysis technique based on 
certain approximations as to how the results are calculated, but with the validity of those 
approximations being checked 

The invention, as embodied below, aims to provide an effective process which 
does not require human expertise to interpret, 

SOURCES OF SPURIOUS INFORMATION 

To understand the operation of the invention it is useful to understand the nature 
of some of the possible error and other information sources which might be accounted for 
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using the present invention. The predominant fonns which can be accounted for are 
contamination with other DNA not originally in the sample; locus or allele drop out for 
information in the sample but not reporting in the results obtained; and stuttering where the 
amplified products include identities which are one repeat unit (or four bases less for a 
teframeric locus) less than the associated allele and may be an allele too or a false 
amplification product. 

Contamination of the DNA sample by DNA which was not originally part of the 
sample is a significant issue when small DNA samples are being amplified, Laboratory 
induced contamination is likely to occur on occasions as equipment and handling is not 
totally clean. Contamination results from fragments of cells in such amplification 
processes would be sufficient to give spurious results. Merely discarding an individual test 
when the negative control produces a result is not a viable option for small DNA sample 
analysis as a proportion of tests might give such a negative control result. 

Locus / allele drop out is also a potential problem particularly with amplification 
of small DNA samples. A heterozygote sample should produce two alleles of the locus 
upon amplification, However, because amplification is an essentially random process, the 
fact that the amplification starts from only a few molecules may mean that problems with 
amplifying one of those alleles at an early stage lead to it not being present to a detectable 
degree in the amplification product This can imply a homozygote identity where in reality 
the identity is heterozygote. 

Stutters are artifacts from short tandem repeat systems and generally represent 
results one repeat unit (or four bases less for a tetiameric locus) than the associated allele, 
Whilst stutters are predictable when large samples are amplified (they generally form a 
15% peak compared with the associated allele peak). This is not the case in small sample 
amplification. As a consequence, stutters can appear close to and even exceed the size of 
the actual allele peaks. This can be a significant issue, particularly if the sample might be 
heterozygote with one of the alleles being four bases less than the other, and could 
consequently be contused with homozygote and stutter result. 



EXPRESSING A LIKELIHOOD OF MATCH BETWEEN TEST SAMPLR AND 
REFERENCE SAMPLES 
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In present analysis systems the result of the expert analysis is either that an 
individual result for the unknown sample is discounted from further consideration or is 
included in the results for further consideration. Thus the raw individual results, set A may 
be whittled down by the expert, excluding individual results which are caused for question, 
to give the individual results included for further consideration, set B in Figure 1 » This 
screened set, set B, is then compared with stored results for various samples of known 
origin, as exemplified by sets C and D in Figure 1 . A match is either agreed between the 
alleles forming the individual results, as between set B and set C in Figure 1, or a match is 
not agreed, as between set B and set D in Figure 1. Thus only some of the initial 
individual results are carried forward and the overall result is a match or not. The extent of 
the question marks over the raw individual results may lead to a substantial number being 
inconclusive and hence the match/non-match decision may be made based on a few points 
only and hence be of reduced statistical significance. 

The technique in this embodiment of the present invention uses a very different 
approach. Once again, Figure 2, the raw results are generated by the analysis process, set 
A. No screening process is carried out, however, and hence no expert input is required. 
Instead all the results are carried forward into set B which is used for comparison purposes 
with the stored results. 

For explanatory purposes only a single stored result is considered for comparison, 
but in reality many such stored results would be considered in an equivalent manner. In 
basic terms the chance that each individual result in the test sample could have arisen 
given the stored sample result for the same locus and the various possibilities by which 
differences between the two could have arisen are considered, Thus, in the case of result E 
from set B alleles 16, 20 are observed and a direct match with the stored sample does not 
occur as that has allele 20 only. However, the technique considers the possibility that the 
16 could have occurred due to contamination of the test sample with DNA have allele 16 
and/or the possibility that the 16 could be a stutter for the 20 allele present in the DNA. In 
a similar vein for result F from set B, allele 22 only, no match occurs as the stored sample 
has alleles 20 and 22. Rather than declaring a non-match, however, the technique 
considers the possibility that the 20 allele was not reported due to allele drop out during 
amplification. 
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In general it is preferred that the result be reported in terms of a likelihood ratio, 
based on the format, likelihood ratio, LR, is given by :- 

LR == Probability of the test sample result if it has the same origin as the stored sample 

Probability of the test sample result if it does not have the same origin as the stored sample 

By pursuing such an approach for each of the individual results of the test sample 
against the individual results of the stored sample an overall likelihood that the test sample 
is a match for that particular stored sample can be obtained. Whilst the LR considers the 
likelihood against any way the test sample result could be reached, this may be restricted 
in practice to reasonable ways to simplify the calculations, but without loss of accuracy as 
the other ways are statistically veTy unlikely. From such overall likelihood ratios for a 
large number of stored samples lists can be obtained of descending likelihood ratios to be 
made between the stored samples and the test sample. Such a listing is more informative 
than a match or non-match statement, is less subjective as it does not involve an initial 
subjective screening process and is potentially more statistically significant as a greater 
number of individual results are considered. 

The above mentioned technique means that all individual results are included for 
consideration; there is no need to exclude feared contamination results, stutter results or 
other artifact results from consideration- These events can be fully accounted for in the 
statistical consideration of the matter. There is also no need to worry that allele drop out 
will lead to an inconclusive or a non-match conclusive as the possibility of this occurring 
can be built in to the consideration. 

Whilst the above mentioned considerations are particularly significant in small 
sample DNA analysis they are potential issues in all such analyses and benefits from the 
use of the technique thus apply in all levels of DNA consideration. The advantages in 
terms of removing subjective decisions also apply whatever the amount of starting DNA. 

This type of technique can significantly simplify the interpretation of mixed DNA 
samples. Such situations occur, for instance, where a suspect may or may not have 
contributed to the mixture along with another known person and/or another unknown 
person as against the sample arising from an unknown person and another known person 
and/or unknown person. As against the previous methodology where interpretation was 
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oflen impossible, the present invention always allows a likelihood ratio to be presented. 
Thus a likelihood ratio may be expressed as > 

LR = Probability of the evidence i f it comes from the stored sample and an unknown 
sample 

Probability of the evidence if it comes from two unknown samples 

In such a case the analysis can consider each stored sample in the database in turn 
in making such an analysis. 

Performing such considerations in an automated manner allows rapid 
consideration of a large number of stored samples against a test sample and/or allows a 
variety of scenarios to be considered (mixture is from suspect A and suspect B; mixture is 
from suspect A and suspect C; mixture is from suspect A and an unknown; etc,) very 
rapidly. This improvement is achieved whilst reducing cost and de-skilling the procedure. 

The manner in which the probability of stored individual results leading to test 
individual result is set out in more detail below, in a different situation, but for which the 
same general principles apply. 

ACCOUNTING FOR ALL RESULTS 

The technique of the invention generally aims to provide a likelihood ratio, LR> 
for the event being considered In many cases this may be the likelihood that a suspect 
vras the source of the DNA being analysed compared to the likelihood that the DNA 
source was someone else. 

Given the potential impact of the other information sources, one or more of these 
is taken into account using the technique of the present invention to provide a fuller 
appreciation of the circumstances which may have generated the DNA sample. In the first 
example discussed the potential impact of contamination and allele dropout axe considered, 
but the technique is capable of application to other information sources and issues too. 

Whatever the likelihood ratio analysis the initial results must be obtained. To 
achieve this the collected sample is subjected to an amplification process such as PCR, to 
make many multiple copies of the DNA present in the initial sample. Where the initial 
sample is formed of only a few cells the number of cycles used may be between 30 and 34 
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cycles to achieve the necessary number of copies. The amplification process is generally 
concerned with producing enough material for effective investigation of the allele identity 
at a large number of locations. Techniques for analysing mixtures are known based around 
the use of short tandem repeats (STR's) as described by Clayton et al (1998) Analysis and 
Interpretation of Mixed Forensic Stains Using DNA STR Profiling, Int. J. Forensic Sci. 91 7 
55-70. 

Example - contamination and allele drop out impact 

In general terms the likelihood ratio in such cases considers the product of :- 

the probability that a given individual result for a locus under test arose, 
given the suspect's identity (homozygote) or identities (heterozygote) for the locus, by the 
various possible routes, including the suspects potential contribution and the potential 
other information source contribution for each of the given individual results; as the 
numerator relative to the product of :- 

the probability that a given individual result for a locus under test arose 
given one of all the possible identity combinations for that locus, by the various possible 
routes, including potential other information source contributions for each of the given 
individual results and each of the possible identity combinations; and 

the probability of that identity combination occurring; as the 

denominator. 

To exemplify this approach in more detail a sample calculation based on three 
different individual results being obtained for a given locus is considered against a 
suspects identity for that locus. The three individual results are :- 

R T -12 
R 2 =16 
R,- 12,16 

with the suspects identity being 12,16. 

If H, is the probability of die evidence if the profile is the suspect's and R 2 is the 
probability of the evidence if the profile is of someone other than the suspect, then the 
process starts by considering all the reasonable identities which might have occurred for 
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the alleles under question. It is possible to account for all the possible identities which 
might have occurred for the alleles, but many of these are so unlikely as to not need 
accounting for in practice. The first column of Table 1 (Figure 3) accounts for the three 
likely identities in this case, 12,12; 12, 16; 16,16. - 

In the next stage the probability of that allele identity occurring is stated, column 
2, giving f i2 2 ; 2f, 2 f, 6 ; f l6 2 ; respectively, the 2f u f 16 reflecting the fact that the identity could 
have been 12,16; or 16,12, This gives the relative balance between the various identities 
being the one to account for one or more of those allele identities being relatively rare, for 
instance. 

In the next stage, column 3, the maimer in which the first individual result, 12, 
could be reached given each of the three allele identities is evaluated. Thus if the allele 
identity where 12, 12, the way in which the individual result 12 occurred would be that 
there was no contamination of the sample and there was no allele drop out for the sample 
in that test. Thus the probability of this individual result arising with this identity is the 
probability that there is no aUele drop out, designated p(D) multiplied by the probability 

that there is no contamination or spurious results, designated p{C) . For this result with 

allele identity 12, 16 the probability is the probability of allele drop out (to account for 16 
not reporting) multiplied by the probability of no allele drop out (to account for 12 
reporting) multiplied by the probability of no spurious reports (to account for no spurious 
16 arising). The equivalent process for result 12 with possible allele identity 16, 16 would 
be the product of the probability of allele dropout (to account for 16 not being reported), 
the probability of spurious reports (to account for 12 being reported) and the probability of 
12 occurring as an allele (to account for differences in the likelihoods of that particular 
spurious allele occurring). 

This process is repeated in columns 4 and 5 of Table 1 for the various individual 
results and possible identity combinations. 

In practice the process would be repeated for a large number of individual results. 

Column 6 of Table 1 represents the product of columns 2, 3, 4 and 5 and relates to 
the overall probability that the allele identity that row represents the set of actual 
individual results. The sum of column 6 gives the overall denominator for the likelihood 
ratio. 
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The numerator for the likelihood ratio is the product of columns 3, 4 and 5 for the 
row having the identity corresponding to the suspect's identity / alleles, in this case row 3. 
Overall this gives the formula 

p(Dfp(cfp{Df 

p{D?p{Qp(D)p(C?fn 2 fK 2 

+ 2p(D) 4 p(C)'p(D) 2 fnf is 

+ p(Dfp{C)p(D)p(C) z fn 2 fu 2 



But this can be reduced, by substitution with p(C), the probability that 
contamination has occurred, and p(C) » the probability that contamination has not 
occurred to :- 



LR 



2fnf\6 



1 + 



fnfi*p(Cf_ 
p(D)p(Df p{Cf 



p(C) can be estimated from a history of observation and/or can be determined for 
a laboratory, for instance, by a series of negative control tests and a consideration of 
spurious bands arising in those. 

p(D) can also be estimated from a history of observation and/or can be determined 
by experimentation. 

Example - evaluation of stutter impact 

As a further example to the accounting for possible contamination and locus/allele 
drop out discussed above this example considers stutter impacts. 
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In the example calculation set out in Table 2 (Figure 4) the position where the 
first individual result gives a,b,c and individual result two gives a,c is considered relative 
to a suspect who is ac and where b is a stutter position. 

Thus in Table 2, column 1 gives the three reasonable possible identities which 
could have given the individual results and column 2 is the probability of those possible 
identities occurring. 

The expressions in column 3 arc derived as follows. Where the result is a 3 b,c and 
the identity being considered is a s b then the result could have occurred p(C) 2 where 

both a and b did not drop out (hence p(D) 2 \ where a spurious report occurs to provide c 
(hence p(C)fJ, and where no stutters for either a or b are reported (hence ). In the next 
row die process generates p(D) 2 (as both a and c report); p(St) (as there is no stutter 

for a); and the combination of p(St).p (C) and p (St) -p(C)^ (to account for the b arising 

as a stutter rather than spurious and the b arising as a spurious not stutter respectively). 

Again the products of columns 2, 3 and 4 for column 5 and the denominator in the 
likelihood ratio is the sum of column 5, the denominator being the product of columns 3 
and 4 of row 3. This gives 

p(Dyp(Qp(Jtf\p(_St)p(Q +P (Jt)p(C)fi] 
2fafifc 2 p(Dfp(D)p(Q 2 p(Sty 

+ Iffotp) <p(C)p(St) 3 [p (St)p(C) + piSt)p(Qfi] 
+ 2fi 2 fiMD) i p(D)p(Q 2 p(Stf 



Again through substitution this can be reduced to > 
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LR « 



1 



2 fife n _p(£) 2 pWfiV° + p( s _!)& 

p(Qp(D){p(St)p(Q + p(St)p(Qfi} 



The types of function used to achieve the likelihood ratio can of course combine 
consideration of contamination, locus / allele dropout and stutter effect in a single case. 
Other such potential sources of information, such as sub-population effects is also possible. 
This is desirable both because it is more consistent with the approaches employed in other 
areas of DNA interpretation but more importantly because it follows logically from the 
correct consideration of the conditional nature of the probability of a genotype of possible 
offenders GIVEN the suspect's genotype. All that is required is to replace p(MV) with the 
conditional probability p^M-IM^ 

Likelihood ratios calculated in this way take into account all of the individual 
results obtained by the analysis process. 

DETERMINING IMPACT OF VARIATIONS IN PROBABILITY OF OTHER 
INFORMATION SOURCES CONTRIBUTING 

As well as enabling the calculation of accurate likelihood ratios in DNA analysis, 
the equations also enable the impact of the likelihood that contamination occurs, dropout 
will occur or stutter will occur upon the analysis process to be considered, Such an 
investigation can be used to determine appropriate thresholds for those probabilities before 
which certain approximations can be deemed to be held true and beyond which certain 
approximations can be deemed to not hold true. This is, potentially pre-ddculated 
information, could be stored and then used in a DNA analysis technique to determine 
whether the results obtained can be processed using an analysing process in which 
approximations for likelihood ratios or other presentation of the results are used. This 
concept is discussed in more detail below together with additional details of the particular 
approximations which might be used in particular circumstances or scenarios- 
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USING APPROXIMATIONS IN ANALYSIS AND RESULTS, WITH CHECK ON VALIDITY 
CONDITIONS FOR THE APPROXIMATION 

Whilst fhe initial embodiment of the invention sets out detailed calculation rates 
for likelihood ratios which take into account each and every one of the results fully, the 
functions used therein can be simplified in certain cases. For instance the formulae 
accounting for potential contamination and potential allele dropout:- 



LR = 



1 



1 + 



fnfiepjC) 2 



p(D)p(Dfp(Q 2 } 



gives a scaling factor of 1 provided p(C) < 0,3 and hence the likelihood ratio approaches 



LR u 



1 

2/12/u 



Thus provided this threshold is not crossed the likelihood ratio can be validly 
calculated based on a substantial simplification. 

Using this process a series of practical situations were considered and the impact 
of contamination and / or locus allele dropout and / or stutter were considered with a view 
to determining the appropriate approximations for such situations and the direshold 
considerations for the applicability of such approximations. 

Example 1A 

Apparently single banded profiles 

When an apparent one-banded homozygotes is encountered in a crime stain 
(Ri=a) and the peak area is small, this may mean that allele dropout has occurred i.e. the 
genotype may in fact be heterozygous: This is considered a possibility whenever the peak 
is close to background. At low peak area, our experimental observation confirms that the 
probability of allele dropout (p(D) is high. If the allele in the crime stain is type a and the 
suspect is type ab then it would seem reasonable to limit Mj to aa, ab or aF wherein F 
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stands for any allele other than type a or type b. This gives, b the general process 
described above:- 



LR = 



1 



2 fa 



2p(D) ' J 



Provided that p(D) > 0.5, the scaling function 



1 



2p(D) J 



> 1,0 as illustrated in 



Figure 5 which is always reasonable when the peak is close to the background, then the 



approximation is demonstrated to be conservative and LR * 



Example IB 

Apparently one banded profiles - effect of one additional replicate 

Taking the previous example, we now consider the advantages of replication, 
where an additional aliquot (R^ of the same DNA extract is separately amplified Suppose 
that the second replicate yields a heterozygote ab that matches the suspect's profile 
(Suspect = ab; R x = a-; R^ab). 

We take account of three possible explanations of the evidence - either M 1 = aa 
homozygote , else M 2 = ab heterozygote. If the first explanation is true then the b allele 
must be a spurious band. In this example the LR would be reported as LR = l/2fa because 
only the a allele was duplicated. The formula that describes this model is: 



LR = 



2fafi 



1 + 



2p(D)p(D)p(Q\ 
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This expression is always less than \l2ff h but the l/2£ evaluation is a 



1 



conservative approximation whenever 



T>1.0 



fi 1 + 



2p(D)p(D)p(C) 



and this is true for all reasonable estimates of p(C) and p(D). 

If a locus appears homozygous with allele % but allele drop-out could have 
occurred, so that the locus was really a heterozygote, then interpretation using l/2£ is 
reasonable provided that contamination is low and the allele peak area itself is small or 
close to the baseline. 

Example 1C 

Additional replicates increase the LR 

Continuing with the previous example, we consider the effect of additional (n) 
replicates that have been analysed and demonstrated to all correspond to the genotype of 
the suspect (ab in this example). The suspect is ab; R 2 ; R 2 _„ = ab (i.e. a total of n+1 
replicates were analysed). The general formula that describes R 3 =a- and n = the number of 
replicates that the genotype ab is: 



Provided that n is greater than or equal to 2, the guideline will allow the reporting 
of LR-V2ff b because both alleles are duplicated. The LR calculated from this equation 
will always be less than this bat is nevertheless a very close approximate (see Figure 6a) 
for most intermediate values of p(D). It is noted that the different between n=2 and u-3 is 
minor in these simulations. Also that the actual values of p(D) and p(C) have very little 
effect on the final estimate (see Figure 6a and 6b) provided that the latter is less than 0.6 
(which it should always be). If n=l then the LR is conservative relative to I/2£. 



LR = 



Ifafh n 



faf^piQ" 



2p(D)[p(D)p(Q\ n 
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Example 2a 

An example where an allegedly contaminant band is observed in conjunction with allele 
dropout 

The next example is more extreme than those previously discussed. Suppose that 
a replicate (RO matches the suspect at one allele (b), but has an additional allele (c) that is 
not found in the suspect under the prosecution hypothesis (H x ). Furthermore, we assume 
that there is no trace of allele a. We assess the condition where the suspect is ab; R x is be; 
R 2 is ab by consideration of the genotypes (Mj) ab, ac, be and bb: 

1 

MD)p(C) _fip(Q 
p(D)p(Q 2p(D)p(D)p(C)_ 




The reporting guideline would only allow reporting of the duplicated b allele 
hence the reported likelihood ratio wouJd be LR = l/2£. There was very little effect 
contributed by p(D) since the scaling function was always greater than 1 ,0 even when p(C) 
was moderately high (see Figure 7a, 7b) demonstrating the conservative nature of the 
reporting guideline. 

Example 2b 

Example where an allegedly contaminant band is observed 

Suppose that the suspect is ab; R } =abc (where c is a supposedly a contaminant 
allele under ,) and R 2 =ab. We limit the possible (Mj) genotypes to ab, ac or be and we 
evaluate against the guideline LR~ V2ff b (see Figures 8a, 8b). Evaluated against p(C)=0.3, 
the approximation is reasonable provided that^ is less than 0.10 and p(D) is less than 0.50. 



Up(D)p(Q 



I p(D)p(Q) 
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Example 3a 

In the case of the stutter counting for formula described above: 



LR = 



1 



2tfc 



1 + 



_ p(e) 2 p(D)fi(/ a + P {st)f c ) 

p(Qp(D){p(St)p(Q + p(St)p(Qfi} 



The scaling function is always less thanl, (as demonstrated by Figure 9) and LR 
approximates to Investigation of this establishes that it is a good approximation 
provided p(St) > 0.5 and p(C) <03. 

Example 3b 

Extension of the dropout definition and evaluation of an actual example 

50FCR replicates of a sample that had 4 homozygous STRs and 6 heterozygous 
loci were analysed. We now expand the definition of p(D) as follows: 

p(D H£ ) - the probability of dropout given that the locus is homozygous. 

p(D H J - the probability that a given allele drops out given that the locos is heterozygous. 

p(D Be2 ) - the probability that both alleles drop out given that the locus is heterozygous. 

In the following calculations the p(D) parameters were either assigned the actual 
observed values in table 3 (Figure 10) or if not available, e.g. (pD H6 ) for locus D8S51 179, 
the mean across available loci was used instead. By observation, p(C) did not exceed 0.2 
for any locus and this value was adopted throughout. Allele frequencies were used from a 
white Caucasian database. 

From two of the replicates, we consider an extreme observation jn table 5 for the 
locus D3S1358, The suspect is 18,18^=15,18 and 1^=15,18. Conventional analysis may 
indicate the results to be either inconclusive or an exclusion since a spurious allele is 
duplicated in the replicates. Using the formula below, limiting our considerations of Mj 
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genotypes to 15,15; 15,18 and 18,18; the LR = 0.068. The evidence supports exclusion, 
but importantly, the LR is greater than zero. 



LR = 



2{p(DHef P (C)} 2 



2+ 



{p{D»o)p(cjf f 15/18 



At locus D8 the suspect is 15 1 5; R1=R2=15 F. The Mj genotypes are limited to 
15 15 and 15 F: 



LR = 



1 



2/is 



/15 {p(DlI e )p(Dse)Y (1-/15) 

2 + P (Dhc) 2 



In HUMTH01, Ri was analysed as 7 F whereas R 2 failed to give a result The 
suspect is 7 9.3. Our evaluation of the LR is limited to 7 F; 7 7 and 7 9.3 genotypes. 
However, the observation that 2 alleles have dropped out in R, is also built into the LR 
calculation below as shown in table 3. 



LR = 



2/7 



1 + 



1 

P (pHo)p{DHo)fl 
1 p{L>He) p(Dffe) p(Dtel) 



Similarly for VWA, Rt felled to give a result and only one allele was observed in 
R>. Our Mj genotype considerations are limited to 19 19 and 19 R The suspect is 19 19: 
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LR = 



2/l9 



1 



fl9 p(Dhc)p(Dh.2)(1- f 19) 

2 p(Dho)p(Dho) 



For Dl 6S539, Rj is 9 12 and " 2 failed to give a result. We limit the evaluation of 
9 9, 9 12 and 12 12 Mj genotypes, considering the possibility of spurious alleles. The 
suspect is 9 12. Again we build the R 2 observation into the LR below; 



LR 



1 



2/9/12 



p(Dm)rtD*)p(&{f*+M 

+ 2p(DHefp(C)p(Dlfe2) 



D18S51 is straight-forward since the 12 16 genotype was observe in both 
replicates. The analysis of D2S133 follows from the equation applied to HUMTH01 
above. 

Examination of the individual LRs calculated using these equations reveals that 
all are either conservative or very close to the estimates derived by calculating LR^2^ for 
homozygotes or LR=2ff b for the D16S539 heterozygote (table 3). When complete locus 
dropout is observed in a replicate this has very little effect, i.e. the scaling function *1.(X 
The combined LR across all loci - 68,000 (using a white Caucasian database) and this 
serves to demonstrate that apparent allele mismatches caused by contamination do not 
necessarily lead to exclusions. 

As is demonstrated by these various examples, therefore, the present invention 
generally provides a technique which can be used to evaluate the impact of variations in 
the probabilities of occurrence of various other information sources and extend that 
information to verification of the accuracy in applying certain assumptions to DNA 
analysis techniques. 



